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Abstract 

This paper establishes Fokker-Planck-Kolmogorov type equations for time-changed 
Gaussian processes. Examples include those equations for a time-changed fractional 
Brownian motion with time-dependent Hurst parameter and for a time-changed Ornstein- 
Uhlenbeck process. The time-change process considered is the inverse of either a stable 
subordinator or a mixture of independent stable subordinators. 

1 Introduction 

A one-dimensional stochastic process X = (X t )t>o is called a Gaussian process if the random 
vector (X tl , • • • , X tm ) has a multivariate Gaussian distribution for all finite sequences < t\ < 
■ ■ ■ < t m < oo. The joint distributions are characterized by the mean function E[A t ] and the 
covariance function Rx{s, t) = Cov(X s , X t ). The class of Gaussian processes contains some of 
the most important stochastic processes in both theoretical and applied probability, including 
Brownian motion, fractional Brownian motion [5, 7, 9, 31], and Volterra processes [1, 8]. 

In this paper, Fokker-Planck-Kolmogorov type equations (FPKEs) for time-changed Gaus- 
sian processes are derived. These FPKEs are partial differential equations (PDEs) satisfied by 
the transition probabilities of those processes. Main features of these equations involve 1) frac- 
tional order derivatives in time and 2) a new class of operators acting on the time variable; see 
(14)-(16). By definition, the fractional order derivative D* of order /? £ (0, 1) in the sense of 
Caputo-Djrbashian is given by 

(i) D ^ = <^ = r^)l!w^ d ^ 

with r(-) being Euler's Gamma function. By convention, set D\ t = d/dt. Introducing the 
fractional integration operator 

J a g(t) = J?g{t) = -L [\t - rr-'g^dr, a > 0, 
r («) Jo 

one can represent t in the form t = ^ o (d/dt) (see [11] for details). 

For the last few decades, time-fractional order FPKEs have appeared as an essential tool for 
the study of dynamics of various complex processes arising in anomalous diffusion in physics 
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[25, 37], finance [13, 17], hydrology [4] and cell biology [30]. Using several different methods, 
many authors derive FPKEs associated with time-changed stochastic processes. For example, 
in [3], FPKEs for time-changed continuous Markov processes arc obtained via the theory of 
semigroups. Paper [14] identifies a wide class of stochastic differential equations whose associ- 
ated FPKEs are represented by time-fractional distributed order pseudo-differential equations. 
The driving processes for these stochastic differential equations are time-changed Levy pro- 
cesses. Two different approaches are taken, one based on the semigroup technique and the 
other on the time-changed Ito formula in [19]. Paper [15] provides FPKEs associated with 
a time-changed fractional Brownian motion from a functional- analytic viewpoint. Continu- 
ous time random walk-based approaches to derivations of time-fractional order FPKEs are 
illustrated in [12, 24, 34, 35]. In the current paper, we follow the method presented in [15]. 

Consider a time-change process E@ = (£"f ) f >o given by the inverse, or the first hitting time 
process, of a /3-stable subordinator W 13 = (wf) t >o with /3 6 (0, 1). The relationship between 
the two processes is expressed as £f = inf{s > 0; W@ > t}. To make precise the problem 
being pursued in this paper, first recall the following results. For proofs of these results as well 
as analysis on the associated classes of stochastic differential equations, see e.g. [12, 14, 15, 24]. 
Throughout the paper, all processes are assumed to start at 0. 

(a) If E@ is independent of an n-dimcnsional Brownian motion B, then B and the time- 
changed Brownian motion {B E p ) have respective transition probabilities p(t, x) and 
q(t,x) satisfying the PDEs 

(2) d t p{t,x) = ^Ap(t,x) and £>£ t q(t, x) = ^Aq(t, x), 

where d t = and A = ^j=i ^ = Sj=i(gfj) 2 ' wrtri the vector x £ 1™ denoted as 



,i 



,x n ). 



(b) Let L be a Levy process whose characteristic function is given by E[e l («> Lf )] = e'^ 
with symbol ip (see [2, 29]). If Er is independent of L, then L and the time-changed 
Levy process (L^s) have respective transition probabilities p(t,x) and q(t,x) satisfying 
the PDEs 

(3) d t p(t,x) = C*p(t,x) and D^ t q(t,x)=C*q(t,x), 
where C* is the conjugate of the pseudo-differential operator with symbol ij). 

(c) If E@ is independent of an n-dimcnsional fractional Brownian motion B H with Hurst 
parameter H £ (0, 1) (definition is given in Example 1), then B H and the time-changed 
fractional Brownian motion {B 1 !^) have respective transition probabilities p(t,x) and 

q(t, x) satisfying the PDEs 

(4) d t p{t, x) = Ht 2H - x Ap(t, x) and D? t q(t, x) = HG p 2H _ lt Aq{t, x), 
where G^t t with 7 G (—1,1) is the operator acting on t given by 



r C+ioo 



(5) G^ t g(t) = m7 + l)Jr P £7±t 



J_ p g(z) 



(*), 



with < C < s and = e^ Ln ^ z \ Ln(z) being the principal value of the complex loga- 
rithmic function ln(z) with cut along the negative real axis. Here g(s) = £t_> s [<?(£)] (s) = 
Ct [g{t)\ (s) and C~\ t [/(s)] (t) denote the Laplace transform and the inverse Laplace trans- 
form, respectively. 
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Obviously (a) provides a special case of both (b) and (c). Note that in (2) and (3) the 
second PDE is obtained upon replacing the first order time derivative dt in the first PDE by 
the fractional order derivative D% t and the right-hand side remains unchanged. On the other 
hand, as (4) shows, this is not the case for a time-changed fractional Brownian motion; the right 
hand side of the time- fractional order FPKE has a different form than that of the FPKE for the 
untime- changed fractional Brownian motion. Namely, the operator G^ H _i t instead of t 211 ^ 1 
appears, which is ascribed to dependence between increments over non-overlapping intervals of 
the fractional Brownian motion B H . When H = 1/2 so that the fractional Brownian motion 
coincides with a usual Brownian motion, the FPKEs in (4) match with those in (2). The 
operators {G^ t ; 7 S (— 1, 1)} are known to satisfy the semigroup property (Proposition 3.6 in 
[15]). 

Following the functional-analytic technique presented in [15] to obtain the time-fractional 
order FPKE in (4), we will establish in Theorem 3 the FPKE for the time-changed Gaus- 
sian process (X E p) under the assumption that the time-change process is independent of 
the Gaussian process X. Moreover, generalization to time-changes which are the inverses of 
mixtures of independent stable subordinators will be considered in Theorem 4. In Section 
4, applications of these results yield FPKEs for time-changed mixed fractional Brownian mo- 
tions as well as those for time-changed Volterra processes. Fractional Brownian motions with 
time-dependent Hurst parameter H = H(t) € (1/2, 1) are among the Volterra processes consid- 
ered; see Example 3. Two equivalent forms of FPKEs for a time-changed Ornstein-Uhlenbcck 
process are compared in Example 5. 

2 Preliminaries 

Let E° be the inverse of a stable subordinator W" starting at with stability index (3 G (0, 1). 
The process W@ is a Levy process with Laplace transform E[e _sM/ * ] = e~ tsf< and self-similarity: 
(Wct)t>o = {c 1 ^ l3 Wf)t>o in the sense of finite dimensional distributions for all c > (see 
[2, 29]). Since W@ is strictly increasing, its inverse E 13 is continuous and nondecreasing, but 
no longer a Levy process (see [22]). The density f E p of i?f can be expressed using the density 
f w? of W? as 

f Bf (r) = d T F(E? <r)=d T {l~ P(Wf < i)} - ~dr{[Jlf w? ] (*)}, 
from which it follows that 

(6) £ t ^ s [f E? (r)}(s) = -d T 

The function f E p (r) is C°° with respect to the two variables t and r. 

The notion of time-change can be extended to the more general case where the time-change 
process is given by the inverse of an arbitrary mixture of independent stable subordinators. Let 
Pi_i(s) = Jq s^dfi(P), where [i is a finite measure with supp/x C (0, 1). Let W 11 be a nonnegative 
stochastic process satisfying E[e~ sW t} = e~ tp ^ and let E$ = inf{r > 0; W>? > t}. Clearly, 
= if (j, = 5p , the Dirac measure on (0, 1) concentrated on a single point /3q. The 
process W^ 1 represents a weighted mixture of independent stable subordinators. Similar to the 
identity (6), the density of has Laplace transform 

(7) C t ^ s [f K (r)](s) = P -Me- T ^ s \ s > 0, r > 0. 



s > 0, r > 0. 
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For further properties of f E f>(r) and /^(r), see [15, 24]. 

The above time-change process is connected with the fractional derivative D* with 
distributed orders given by 

(8) DSg(t) = D^ t g(t) = [ DPg(t)d»(p). 

Jo 

Namely, if is replaced by E^ in items (a) and (b) of the list of known results in Section 
1, then the FPKEs for the time-changed Brownian motion (B E ») and the time-changed Levy 
process {L E v) are respectively given by (see [14]) 

Dt, t q{t,x) = ^Aq(t,x) and D^ t q(t,x) = C*q(t,x). 

For investigations into PDEs with fractional derivatives with distributed orders, see [20, 23, 33]. 

The covariancc function of a given zero-mean Gaussian process is symmetric and positive 
semi-definite; conversely, every symmetric, positive semi-definite function on [0, oo) x [0, oo) 
is the covariance function of some zero-mean Gaussian process (see e.g., Theorem 8.2 of [18]). 
Examples of such functions include Rx(s, t) = sAt for Brownian motion and Rx(s, t) = a$-\-s-t 
which is obtained from linear regression (see [28]). The sum and the product of two covariance 
functions for Gaussian processes are again covariance functions for some Gaussian processes. 
For more examples of covariance functions, consult e.g. [28]. 

An n-dimensional Gaussian process X = (X 1 , . . . , X n ) is a process whose components X 3 
are independent one-dimensional Gaussian processes with possibly distinct covariance func- 
tions Rxj(s,t). The variance functions R X i{t) — Rxi(t,t) will play an important role in 
establishing FPKEs for Gaussian and time-changed Gaussian processes. For differentiablc 
variance functions R X i(t), let 

1 " 

(9) A = A x = -Y,R'xMd%- 

3 = 1 

Clearly A = ^ A if X is an n-dimensional Brownian motion. 

Proposition 1. Let X = [X , ...,X n ) be an n-dimensional zero-mean Gaussian process 
with covariance functions Rxi(s,t), j = l,...,n. Suppose the variance functions Rxi(t) = 
Rx^(t,t) are differentiablc on (0, oo). Then the transition probabilities p(t,x) of X satisfy the 
PDE 



(10) d t p(t, x) = Ap(t, x), t>0, 

where A is the operator in (9). 

Proof. Since the components X 3 of X are independent zero-mean Gaussian processes, it follows 
that 



P (t, x ) = n(2^(*)r i/2 x cx P 



Direct computation of partial derivatives of p(t,x) yields the equality in (10). □ 
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Remark 2. a) In Proposition 1, if the components X 7 are independent Gaussian processes 
with a common variance function Rx{t) which is differentiable, then (10) reduces to the 
following form which agrees with the classical FPKE: 

(12) d t p(t,x) = ^R' x (t)Ap(t,x). 

b) If X = (X 1 , . . . , X n ) is an n-dimcnsional Gaussian process with mean functions m X j (t) 
and covariancc functions Rxj(s,t), and if both mxiif) and Rxj(t) = Rxi{t,t) are differen- 
tiable, then the associated FPKE contains an additional term: 

n 

(13) 9 t p(t, x) = Ap(t, x) + Bp(t, x) where B = — Y"^ m' xj (t) d x j . 

i=i 

Such Gaussian processes include e.g. the process defined by the sum of a Brownian motion and 
a deterministic differentiable function. 

c) The initial distribution of the Gaussian process X in Proposition 1 is not given; however, 
it needs to be specified in order to guarantee uniqueness of the solution to PDE (10). The 
same argument applies to the PDEs to be established in Theorems 3 and 4 as well. 



3 FPKEs for time-changed Gaussian processes 

Theorems 3 and 4 formulate the FPKE for a time-changed Gaussian process under the as- 
sumption that the time-change process is independent of the Gaussian process. The case 
where those processes are dependent is not discussed in this paper. As in Proposition 1, the 
variance function plays a key role here. 

Theorem 3. Let X = (X ,...,X n ) be an n- dimensional zero-mean Gaussian process with 
covariancc functions Rxi (s, t), j = 1, . . . , n, and let E" be the inverse of a stable subordinator 
of index (3 £ (0,1), independent of X . Suppose the variance functions Rxi{t) = Rxi{t,t) 
are differentiable on (0, oo) and Laplace transformable. Then the transition probabilities q(t, x) 
of the time-changed Gaussian process (X E /i) satisfy the equivalent PDEs 

n 

(14) D^ t q(t,x) ^Y. J t~ Pk x^t d l^ x ^ * > °> x e R ™ 

3=1 
n 

(15) d t q(t,x) =J2A p xi)t e^q{t,x), t > 0, x € K", 

3=1 

where A xj j = 1, . . . , n, are the operators acting on t given by 



(16) A X3t9 (t) = ^C^ t 



/ ^-z^)R xi {^-zP)g{z)dz 

2m 



(t), 



with z 13 — e^ Ln ( z ), Ln(z) being the principal value of the complex logarithmic function ln(z) 
with cut along the negative real axis, and C being a curve in the^complex plane obtained via 
the transformation £ = z@ which leaves all the singularities of Rxi on one side. 

Proof. Let p(t,x) denote the transition probabilities of the Gaussian process X. For each 
x £ K™, it follows from the independence assumption between E" and X that 

/•OO 

(17) q(t,x)= f E f,(T)p{T,x)dT, t>0. 
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Relationship (17) and equality (6) together yield 

(18) q{s, x) = s^pis 13 , x), s>0. 

Since Rxo (t) is Laplace transformable, Rxj (s) exists for all s > a, for some constant a > 0. 
Taking Laplace transforms on both sides of (10), 



1 ™ 

(19) sp(s,x) -p(0,x) = -J2^s[R' x j(t)d%p(t,x)}(s) 

3=1 
1 ™ 

= -^(A[i?^(t)]*£ t [e-K*^)])( s ) 



3=1 



2 ^ 

3=1 

2 ^ 

3=1 



i />c-\-ioc 

— R> XJ (s-0d 2 x jP(Cx)dC 



l 



— / Bf xj (8-z")a^p(zP,x)z 



s > a, 



where * denotes the convolution of Laplace images and the function 
(20) R^ j (s) = sR x ' j (s), s>a, 

exists by assumption. Equation (20) is valid since R X j (0) = due to the initial condition 
X?(0) = 0. Since Eq = with probability one, it follows that p(0,x) — q(0,x). Replacing s 
by s 13 and using the identity (18) yields 

P ' 



(21) 



sq(s,x) - q(0,x) 



E 

3 = 1 



1 

2TTi 



R'{s -zP)d 2 x3 q(z,x)dz 



, s > a 



1//3 



Since the left hand side equals Ct^ s [dtQ(t,x)](s), PDE (15) follows upon substituting (20) and 
taking the inverse Laplace transform on both sides. Moreover, applying the fractional integral 
operator j\~ to both sides of (15) yields (14). □ 

The next theorem extends the previous theorem to time-changes which are the inverses of 
mixtures of independent stable subordinators. 

Theorem 4. Let X = {X 1 ,...,X n ) be an n- dimensional zero-mean Gaussian process with 
covariance functions Rxi{s,t), j = l,...,n, and let E^ be the inverse of a mixture of 
independent stable subordinators, independent of X . Suppose the variance functions Rxi{t) = 
Rxi{t,t) are differ entiable on (0,oo) and Laplace transformable. Then the transition probabil- 
ities q(t,x) of the time-changed Gaussian process (X E />) satisfy the PDEs 



(22) 

and 
(23) 



n -i 

D^ t q(t,x) = V / Jl- p K X3 t d 2 xJ q(t,x)dfi(f3), t > 0, x G 
3=1 J ° 



d t q(t,x) = J2Ki,t d xil( t ' x )> <>0 > xeR n , 

3=1 
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where D% t is the operator in (8) and A xj j = 1, . 



, n, are the operators acting on t given 



by 



(24) A^.g(t) = -C~ s 



^— / Ov( s ) ~ Pn( z )) R xi{Pn{s) - p ll {z))m ll {z)g{z)dz 



(*), 



iwif/i Pfi(z) = Jq 1 e' 31 " 11 '^ dp((3), m^z) = - /3z@ dp((3) : and C being a curve in the complex 

plane obtained via the transformation £ = pft(z) which leaves all the singularities of Rxi on 
one side. 

Proof. We only sketch the proof since it is similar to the proof of Theorem 3. Let p(t,x) 
denote the transition probabilities of the Gaussian process X. For each x £ K™, it follows from 
relationship (17) with f E $ replaced by f E p> together with equality (7) that 



(25) 



s > 0. 



Taking Laplace transforms on both sides of (10) leads to the second to last equality in (19). 
Letting ( = p^(z) yields 



sp(s,x) -p(0,x) = - 



3=1 



1 

2ttI 



p fl {z))d 2 x] p(p fl (z),x)^^-m fl {z)dz 



which is valid for all s for which Rxj(s) exists. Replacing s by Pa(s) and using the identity 
(25) yields an equation similar to (21). PDE (23) is obtained upon taking the inverse Laplace 
transform on both sides. Finally, applying the fractional integral operator J t 1_/3 and integrating 
with respect to p, on both sides of (23) yields (22). □ 

Remark 5. a) If p = 5, 3a with f3 e (0, 1), then A£ 3 t g(t) = A^ t g(t) and the FPKEs in (22) 
and (23) respectively reduce to the FPKEs in (14) and (15) with (3 = (3q, as expected. 

b) In Theorem 4, if the components X J are independent Gaussian processes with a common 
variance function Rx[t) which is diffcrentiable and Laplace transformable, then the FPKEs in 
(22) and (23) respectively reduce to the following simple forms: 



(26) 
(27) 



4- p A xt Aq(t,x)dp(/3); 



d t q(t,x)=A xt Aq(t,x) 



c) As Example 1 in Section 4 shows, Theorem 3 extends the time-fractional order FPKE in 
(4) for a time-changed fractional Brownian motion to that for a general time-changed Gaussian 
process, revealing the role of the variance function in describing the dynamics of the process. 



4 Applications 

This section is devoted to applications of the results established in this paper concerning 
FPKEs for Gaussian and time-changed Gaussian processes. For simplicity of discussion, we 
will consider the time-changed process given by the inverse of a stable subordinator W" 
of index f3 £ (0, 1), rather than the more general time-change process . 

Example 1. Fractional Brownian motion. One of the most important Gaussian processes in 
applied probability is a fractional Brownian motion B H with Hurst parameter H £ (0, 1). A 
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one-dimensional fractional Brownian motion is a zero-mean Gaussian process with covariance 
function 

(28) R BB (s,t)=E[B?B*] = \{ S 2H +t 2H - \s-t\ 2H ). 

If H = 1/2, then B H becomes a usual Brownian motion. An n-dimcnsional fractional Brownian 
motion is an n-dimensional process whose components arc independent fractional Brownian 
motions with a common Hurst parameter. 

Stochastic processes driven by a fractional Brownian motion B H are of increasing interest 
for both theorists and applied researchers due to their wide application in fields such as math- 
ematical finance [7], solar activities [31] and turbulence [9]. The process B , like the usual 
Brownian motion, has nowhere diffcrcntiablc paths and stationary increments; however, it docs 
not have independent increments. B H has the integral representation Bf = j Q Ku{t, s)dB s , 
where B is a Brownian motion and Kjj(t, s) is a deterministic kernel. B H is not a semimartin- 
gale unless H = 1/2, so the usual Ito's stochastic calculus is not valid. For details of the above 
properties, see [5, 27]. 

Let B H be an n-dimensional fractional Brownian motion and let E° be the inverse of a 
stable subordinator of index (3 £ (0,1), independent of B H . Then the components of B H 
share the common variance function (t) = i 2H and its Laplace transform R' BH ( s ) = 
2HT(2H)/s 2H . Hence, Proposition 1 and Theorem 3 immediately recover both the FPKEs in 
(4) for the fractional Brownian motion B H and the time-changed fractional Brownian motion 
(B^/j). In this case, 

(29) j\ ^ B H it = HG^ H _ 1V 

where t is the operator given in (5). Note that the curve C appearing in the expression of 

the operator Ag H , in (16) can be replaced by a vertical line {C + ir; r el} with < C < s 
since the integrand has a singularity only at z = s. 

Example 2. Mixed fractional Brownian motion. Let X be an n-dimcnsional process de- 
fined by a finite linear combination of independent zero- mean Gaussian processes X\, . . . , X m : 
Xt = YliLi a e^e,t with a\, . . . , a rn £ R. For simplicity, assume that for each I = 1, . . . , m, 
the components of the vector Xg = (Xl, . . . , X™) share a common variance function Rx e (t). 
Then the process X is again a Gaussian process whose components have the same variance 
function Rx(t) = Sfli a ?-^Q Therefore, it follows from Proposition 1 and Theorem 3 
that the FPKEs for X and (X E f>), under the independence assumption between E@ and X, 
are respectively given by 

(30) dtp(t,x) = <f>(t)Ap(t,x) and D^ t q{t,x) = $f Aq(t,x), 

where 4>{t) = | J2T=i a l^'x e (^) an d ®t = XX=i a lJt~^ A Xt t . Notice that <fi(t) simply denotes 
the multiplication by a function oit whereas $f is an operator acting on t. This generalizes the 
correspondence between the function t 2H ~ x and the operator G^ H _ X t observed in the FPKEs 
in (4) for the fractional Brownian motion and the time-changed fractional Brownian motion. 

A mixed fractional Brownian motion is a finite linear combination of independent fractional 
Brownian motions (see [26, 32] for its properties). It was introduced in [6] to discuss the price 
of a European call option on an asset driven by the process. The process X considered in that 
paper is of the form X t = B t + a Bp , where a £ R, B is a Brownian motion, and B H is a 
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fractional Brownian motion with Hurst parameter H £ (0, 1). In this situation, the FPKEs in 

(30) , with the help of (29), yield 

(31) d tP (t, x) = - Ap(t, x) + a 2 Hi 2 "- 1 Ap(t, x); 

(32) D^ t q{t,x) = ^Aq(t,x) +a 2 HG p 2H _ ht Aq(t,x). 

Example 3. Fractional Brownian motion with variable Hurst parameter. Volterra processes 
form an important subclass of Gaussian processes. They are continuous zero-mean Gaussian 
processes V = (Vt) defined on a given finite interval [0,T] with integral representations of the 
form V t = y K(t,s)dB s for some deterministic kernel K (t, s) and Brownian motion B (see 
[1, 8] for details). Fractional Brownian motions are clearly an example of a Volterra process. 
In particular, if B H is a fractional Brownian motion with Hurst parameter H £ (1/2, 1), then 
B H = J*K H (t,s)dB s with the kernel 

(33) K H {t,s) = c H s 1 ' 2 - H j\r-s) H -^ 2 r H - x f 2 dr, t>s, 

where the positive constant ch is chosen so that the integral J* As Kn(t, t)Kh (s, r) dr coincides 
with R B H(s,t) in (28). Increments of B H exhibit long range dependence. 

A particular interesting Volterra process is the fractional Brownian motion with time- 
dependent Hurst parameter H(t) suggested in Theorem 9 of [8]. Namely, suppose H{t) : 
[0,T] — > (1/2, 1) is a deterministic function satisfying the following conditions: 

(34) inf H(t) > \ and H(t) £ S 1/2+a ,2 for some a £ ( 0, inf H(t) - - 

where S v ,2 is the Sobolev-Slobodetzki space given by the closure of the space C 1 [0, T] with 
respect to the semi-norm 

Then representation (33) with H replaced by H{t) induces a covariance function Ry(s,t) = 
Jq* As Kfj(t){t,r)Kf{( s )(s,r)dr for some Volterra process V on [0,T]. The variance function is 
given by Rv(t) = t 2H ^ and is necessarily continuous due to the Sobolev embedding theorem, 
which says S Vi 2 C C[0,T] for all r\ > 1/2 (see e.g. [16]). Therefore, H(t) is also continuous. 

Let H(t) : [0, oo) — > (1/2,1) be a differentiable function whose restriction to any finite 
interval [0,T] satisfies the conditions in (34). For each T > 0, let K V T(s,t) be the kernel 
inducing the covariance function R V T(s,t) of the associated Volterra process V T defined on 
[0, T] as above. The definition of K V T(s,t) is consistent; i.e. K V T 1 (s,t) = K V T 2 (s,t) for any 
< s, t < Ti < T2 < 00. Hence, so is that of R v t (s, t), which implies that the function Rx (s, t) 
given by Rx(s, t) = Ryr (s, t) whenever < s, t < T < 00 is a well-defined covariance function 
of a Gaussian process X on [0, 00) whose restriction to each interval [0,T] coincides with 
V T . The process X represents a fractional Brownian motion with variable Hurst parameter. 
The variance function Rxif) = t 2H ^ is differentiable on (0, 00) by assumption and Laplace 
transformable due to the estimate Rx (t) < t 2 . Therefore, Proposition 1 and Theorem 3 can be 
applied to yield the FPKEs for X and the time-changed process (X E p ) under the independence 

assumption between E@ and X. 
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Example 4. Fractional Brownian motion with piecewise constant Hurst parameter. The 
fractional Brownian motion discussed in Example 3 has a continuously varying Hurst pa- 
rameter H(t) : [0,oo) — > (1/2,1). Here we consider a piecewise constant Hurst parameter 
H(t) : [0, 00) — > (0, 1) which is described as 

N 

(36) H(t)=J2HkI [Tk ,T k+1 M 

where {H k }^ =0 are constants in (0, 1), {T k }^ =0 are fixed times such that = To < 7\ < • • • < 
Xjv < Tn+i = 00, and I[T k ,T k+1 ) denotes the indicator function over the interval [Tfc,Tfc + i). 

For each k = 0, . . . , N, let B Hk be an n-dimcnsional fractional Brownian motion with Hurst 
parameter H k . Let X be the process defined by 

fe-i 

(37) X t = ( B ?- +1 - B ?- ) + i, B t" - B T k ) whenever t G [T k , T k+1 ). 

j=o 

Then X is a continuous process representing a fractional Brownian motion which involves 
finitely many changes of mode of Hurst parameter (described in (36)). 

The transition probabilities of the process X are constructed as follows. For each k = 
0,...,N, let 6 k (t) = Hkt 211 "- 1 for t G [T k ,T k+1 ). Let {p k (t, x)}f =0 be a sequence of the 
unique solutions to the following initial value problems, each defined on [T k ,T k +i) x M n : 

(38) d tPo {t,x) = e (t)A Po (t,x), tG(0,Ti), 

(39) p {0,x) = 6o{x), x£R n , 

where Sq(x) is the Dirac delta function with mass on 0, and for k = 1, . . . , N, 

(40) d tPk (t,x) = k (t)A Pk {t,x) 7 t e (T k ,T k+1 ), x G K n , 

(41) p k (T k> x)=p k - 1 {T k -0,x), x G R". 

Define functions 8(t) and p(t,x) respectively by 9(t) = 9 k (t) and p(t,x) = p k (t,x) whenever 
t G [T k ,T k+ i). Then the transition probabilities of X axe given by p(t,x) and satisfy 

(42) d t p(t,x) = 9(t)A P (t,x), t G UL)(7fc> T fc+i)> * e M n , 

(43) p(0,x) =<5 (x), a; G E™. 

(44) p(T fc) a;) =p(T fc -0,x), k = l,...,N. 

Discussion on existence and uniqueness of the solution to this type of initial value problems is 
found in [36] . To ensure transition probabilities which are continuous in time, the initial value 
problem associated with the time-changed process (X E p ) is given by 

(45) D P „ t q{t,x) = <d^Aq(t,x), t G (0,oo), x G M", 

(46) q{0,x) = 5 (x), x G R n , 

(47) q(T k , x) = q(T k — 0, x), x G I", k = 1, . . . ,N, 

where Gf is the operator acting on t defined by 0f = X^fclo ^fc^2_f/ fc -i t^[T k ,T k+1 ){t)- 

Remark 6. Combining ideas in Examples 3 and 4, it is possible to construct a fractional 
Brownian motion having variable Hurst parameter H(t) G (1/2, 1) with finitely many changes 
of mode and to establish the associated FPKEs. 
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Example 5. Ornstein-Uhlenbeck process. Consider the one-dimensional Ornstein-Uhlenbeck 
process Y given by 

(48) Y t = y e- at + a f e - a ^~ s UB s , f > 0, 

Jo 

where a > 0, a > 0, yo G M. arc constants and B is a standard Brownian motion. If a = 0, 
then Yt = yo + uBt, a Brownian motion multiplied by a starting at yo. Suppose a > 0. The 
process Y defined by (48) is the unique strong solution to the inhomogeneous linear SDE 

(49) dY t = -aY t dt + adB t with Yo = yo, 
which is associated with the SDE 

(50) dY t = -aY t dE^ + adB E ,3 with Y Q = y , 

via the dual relationships Y t = Y E p and Y t = Y w n ; see [19] for details. 
Consider the zero-mean process X defined by 

(51) X t = Y t -y e- at = a [ e^'" 8 ) dB s . 



X is a Gaussian process since each random variable X t is a linear transformation of the 
Ito stochastic integral of the deterministic integrand e as . Direct calculation yields Rx(t) = 
§jj(l — e~ 2at ) and R' x ( s ) — s +2a • Therefore, due to Proposition 1 and Theorem 3, the initial 
value problems associated with X and (X E n), where E@ is independent of X, are respectively 
given by 

2 

(52) d t p(t,x) = ye- 2ai ^p(i,4 P(0,x) = 5 (x); 



s-tt 



(53) D ^ t q(t,x) = ^J t L - p £ 



dz 



(t), q(0,x) = 6 (x). 



The unique representation of the solution to the initial value problem (52) is obtained via the 
usual technique using the Fourier transform. Moreover, expression (17) guarantees uniqueness 
of the solution to (53) as well. 

Notice that the two processes X and {X E n ) are unique strong solutions to SDEs (49) and 
(50) with yo = 0, respectively. Therefore, it is also possible to apply Theorem 4.1 of [14] 
to obtain the following forms of initial value problems which are understood in the sense of 
generalized functions: 

(54) d t p(t,x) = ad x {xp(t,x)} + —d*p(t,x), p(0, x) = 5 (x); 

2 

(55) D^ t q{t,x) = ad x {xq{t,x)} + y ^g(t,s), q(0,x) = 5 (x). 

Actually these FPKEs hold in the strong sense as well. For uniqueness of solutions to (54) and 
(55), see e.g. [10] and Corollary 3.2 of [14]. 

The above discussion yields the following two sets of equivalent initial value problems: (52) 
and (54), and (53) and (55). At first glance, PDE (52) might seem simpler or computationally 
more tractable than PDE (54); however, PDE (53) which is associated with the time-changed 
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process has a more complicated form than PDE (55). A significant difference between PDEs 
(52) and (54) is the fact that the right-hand side of (54) can be expressed as A*p(t, x) with 
the spatial operator A* = ad x x + \d^. whereas the right-hand side of (52) involves both the 
spatial operator d% and the time- dependent multiplication operator by e~ 2at . This observation 
suggests: 1) establishing FPKEs for time-changed processes via several different forms of 
FPKEs for the corresponding untime-changed processes, and 2) choosing appropriate forms 
for handling specific problems. 

Remark 7. Example 5 treated a Gaussian process given by the ltd integral of the deterministic 
integrand e~ atyt ~ s ^ . Malliavin calculus, which is valid for an arbitrary Gaussian integrator, can 
be regarded as an extension of Ito integration (see [18, 21, 27]). It is known that Malliavin- 
type stochastic integrals of deterministic integrands are again Gaussian processes. Therefore, 
if the variance function of such a stochastic integral satisfies the technical conditions specified 
in Theorem 3, then the FPKE for the time-changed stochastic integral is explicitly given by 
(14), or equivalently, (15). 
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